Photon scattering from strongly driven atomic ensembles 
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The second order correlation function for ligiit emitted from a strongly and near-resonantly driven 
dilute cloud of atoms is discussed. Because of the strong driving, the fluorescence spectrum separates 
into distinct peaks, for which the spectral properties can be defined individually. It is shown that the 
second-order correlations for various combinations of photons from different spectral lines exhibit 
bunching together with super- or sub-Poissonian photon statistics, tunable by the choice of the 
detector positions. Additionally, a Cauchy-Schwarz inequality is violated for photons emitted from 
particular spectral bands. The emitted light intensity is proportional to the square of the number 
of particles, and thus can potentially be intense. Three different averaging procedures to model 
ensemble disorder are compared. 

PACS numbers: 42.25.Hz, 42.50.Hz, 42.50.Ar, 42.50.Dv 



I. INTRODUCTION 

The characterization of light is an ubiquitous problem, 
and a convenient formalism for this is the use of correla- 
tion functions. However it is well known that it is not pos- 
sible to sufficiently distinguish the nature of a given light 
source from only the first-order correlation function 111]. 
In particular, the quantum properties of light cannot be 
extracted from the first order correlation function. This 
motivated the study of second-order correlations, initi- 
ated by the intensity-correlation experiments conducted 
by Hanbury-Brown and Twiss . Subsequently, second- 
order correlation measurements have found applications 
in many fields of modern physics [3| such as astronomy 
optics 1 5], high-energy physics jg], condensed matter 
physics [7, ^] and atomic physics . 

Next to the characterization of a physical system via 
higher-order correlations, it has also been suggested to 
use strongly correlated particles or even entangled par- 
ticles as input to a system for various applications, with 
the most obvious examples of quantum computation, 
quantum communication and quantum information pro- 
cessing. This prompts the question for efficient sources 
of correlated photons. A standard method for generat- 
ing entangled photons is the parametric down-conversion 
processes [10|. Alternatively, entangled photons can be 
generated in four- wave mixin g [ 111 Il2| or electromagnet- 
ically induced transparency [ij] processes. Further, an 
experiment on photon anti-bunching in phase-matched 
multi-atom resonance fluorescence was reported in [l^ . 
Non-classical photon pairs for scalable quantum commu- 
nication with atomic ensembles and ultraviolet entangled 
photons in a semiconductor were generated in [15] and 
[la ] , respectively. An atomic memory for correlated pho- 



ton states was demonstrated in llTl. and a heralded en- 
tanglement source was realized in [I^. However, despite 
the large variety of sources, in many cases applications 
are restricted by a limited production rate for correlated 
photons, and by the lack of sources with tunable corre- 
lations. 

Interestingly, it has been shown that disordered me- 
dia can sometimes exhibit surprising correlation or coher- 
ence effects. For example, speckle patterns in the quan- 
tum correlations within li ght scattered by a disordered 
medium were observed in [l^. A cooperative radiation 
force in the presence of disorder was observed as well [l^l . 
An experimental investigation of the origin of disorder in 
parametrically excited waves on a fluid surface (Fara- 
day waves) was performed in [2l| . Instabilities of waves 
in nonlinear disordered media and subwavelength spatial 
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FIG. 1: (Color online) Schematic setup of a dilute atomic 
ensemble pumped with a coherent field with wave- vector fe_L. 
(a) shows the energy levels of each of the ensemble particles, 
and the interaction with the strong coherent light with cou- 
pling strength fi. (b) depicts the ensemble with typical inter- 
particle distance length scale denoted by I with I ^ 2n/kL- 
We consider the case of photon pair emission in forward direc- 
tion and denote the angle between the two emitted photons 
with wave- vectors k\ and as (^o, and fci -|- ^2 ~ 2^^. The 
direction of the emission cone defined by fci and is charac- 
terized by the angle (j) between fci -|- and kh ■ 
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correlations in near-field speckle patterns were discussed 
in [2^ [2^. Correlation and coherence can also be ob- 
served in the context of localization Furthermore, 
coherent backscattering was investigated in [25| , and ran- 
dom lasing in [26. - ,28,] . 

In this paper, we study the generation of correlated 
light from a disordered many-particle system, and in par- 
ticular focus on the case of a strong near-resonant driving 
of the ensemble. Due to the strong driving, the fluores- 
cence spectrum dissolves into distinct lines, which can 
be characterized separately. We then investigate second- 
order correlation functions of light scattered into the var- 
ious spectral lines. The disordered ensemble is modelled 
by calculating the emitted light for a fixed orientation 
between the atoms, followed by a suitable configuration 
averaging. For this, we compare three different averag- 
ing procedures. We find that the second-order correla- 
tion functions typically exhibit strong correlations if de- 
tected close to the forward scattering direction, whereas 
the corresponding first-order correlation functions are 
isotropic. In particular, bunching together with super- or 
sub-Poissonian photon statistics can be achieved for dif- 
ferent combinations of spectral lines, and in some cases 
the statistics is tunable by the choice of the detector po- 
sitions. We additionally demonstrate that light emitted 
from certain spectral bands violates a Cauchy-Schwarz 
inequality. Since the unnormalized second-order correla- 
tion functions of the emitted light scales with the number 
of particle in the ensemble squared, intense beams of cor- 
related photons can be generated. 



II. THEORETICAL CONSIDERATIONS 
A. The model 

We proceed by investigating an atomic sample of ar- 
bitrary shape and of characteristic size d, consisting of 
distinguishable non-overlapping two-level particles. We 
assume conditions such that multiple scattering effects 
can be neglected, i.e., we consider the single scattering 
regime in which the photon mean free path length is of 
the order of the sample size or larger. Related, we also 
restrict the typical inter-particle separation I to satisfy 

AL<c;<d, (1) 

with d/c < Tg, where is the spontaneous decay time. 
All particles have identical atomic transition frequen- 
cies ujQ, and are localized at random positions rj with 
j G {1,2,- •• ,N}. We define the interparticle separa- 
tion vectors as rij = fi — fj. The external laser field 

has frequency ujl = ck^ = 2t:c/Xl, wave vector fc^ and 
wavelength Al (see Fig. [1]). 

Under the action of the laser field, the system is best 
described in a suitable dressed state picture. In elec- 
tric dipole and rotating wave approximations, the system 



Hamiltonian can be written a.s H = Hq + Hi , where [2^ 

N 

Hq =^h{ujk - i^L)alak +'^fiCljR^j , (2a) 

N 

HI^^Y.Y.(9k■dM4S-e-'^^'-~'-^■'^~ll.c.}, (2b) 

k j=i 

S- = ^ sin 20 J - Ri{^ sin^ 0^ + i?^^' cos^ 0^ . (2c) 

Here, Hq represents the Hamiltonian of the free electro- 
magnetic field (EMF) and free dressed atomic subsys- 
tems, respectively, while Hj accounts for the interaction 
of the laser-dressed atoms with the EMF. and are 
the field annihilation and creation operators obeying the 
standard commutation relations for bosons. The atomic 
operators R^^^ = \a)jj{P\ describe the transitions be- 
tween the dressed states |/3)j and \a)j in atom j for a ^ /3 
and dressed-state populations for a = (3, and satisfy the 
commutation relations of the su(2) algebra. The dressed 
states \a)j entering the operators R^^^ can be represented 
through the bare states \a)j via the transformations 

|l)j = sin6l|2)j +cos6l|i)j , (3a) 
|2)j =cos6'|2)j -sin6l|i)j. (3b) 

We further defined 

i?., = |2),,(2|-|i)„-(i|, (4) 

which is the difference of the upper and lower dressed 
state population. Further, 

n^nj = yW+JK/W (5) 

is the generalized Rabi frequency, with 2n — {d ■ EL)/h. 
Here, El is the electric laser field strength, and d = dj 
is the transition dipole matrix element. The detuning 
Is. = uq ~ UJL is characterized by cot 20 ~ A/(2r2). 

B. Spectral decomposition 

In the following we make use of the fact that the reso- 
nance fluorescence spectrum of the light scattered by an 
atomic system pumped by a strong near-resonant driving 
field splits up into distinct lines, which are known as the 
Mollow spectrum [s^l- The condition for well- resolved 
spectral lines is J7 » 7, with 7 = 1/ts being the single- 
atom spontaneous decay rate. In this limit of strong driv- 
ing, it is reasonable to define optical properties for each 
of the spectral lines separately. Thus, it follows from the 
interaction Hamiltonian Eq. (j2bl) that the different oper- 
ators in Eq. (Pel) can be considered as sources of individ- 
ual spectral lines. In particular, in the following we will 
denote light originating from Rzj sin(20)/2 as the central 

spectral component indicated by C, and R21 cos^ and 
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r\^2 sii^^ ^ the right (R) and left (L) spectral sideband 
components. These sidebands are emitted at frequencies 
u!^ = uj[^ + 2^1 and a;_ = — 2fi, respectively [sTI. [s^. 
In what follows, we shall use this decomposition to inves- 
tigate the properties of the scattered light. 



C. First order correlations 

In order to calculate the light scattered by the cloud of 
two-level scatterers, we assume detection in the far-zone 
limit, that is, the linear dimension of the atomic system 
d is much smaller than the distances between the cloud's 
center-of-mass and the detector at R. The intensity of 
the scattered light can be calculated from the first order 
photon correlation function as 



I„^iR) = {al{R)a,r^iR)) 



(6) 



where photon operator for the mth spectral band 

with m 6 {C, R, L}. If was shown in [23] that these 
intensities of the different spectral lines can be evaluated 



Ic{R) '^RirJ^,LUL) (R.jR.^) sin2(20) , (7a) 

N 

Il{R)^J2 (^^'-^ ' ) (^12 4i ) sin^ e, (7b) 

N 

Here, 'ifii{rji,Lu) — 'i' ii{uj) exp[i{k ~ kL)rji] for a scat- 
tered photon of wave- vector k. 5'/j(a;) depends on the 
atom-environment coupling and in general is a function 
of frequency with R = \R\ :§> k^^. As we focus on large 
inter-particle separations [see Eq. ([T])] , there is only a neg- 
ligible direct coupling between the atoms, and thus the 
collective correlators {R12R21) describing the vacuum- 
mediated interactions among the emitters i and j decou- 
ple for i ^ that is. 



\R12R21 ) 



(Rihi^i^i) 



(* 7^ j) ■ (8) 



In the strong driving field case, the atomic variable {R'f']) 
scales as j/d and therefore can be neglected in a secu- 
lar approximation. As a consequence, the intensities of 
the spectral side-bands are proportional to the number 
of scatters in the ensemble, i.e., {Il,Ir} oc N. Near the 
resonance {9 = 7r/4), the intensity of the central spectral 
line is also proportional to TV. In what follows, we will 
restrict the analysis to the resonant driving field case in 
which the intensities of all spectral lines are proportional 
to the number of scatters N. In this case, the intensities 
of spontaneously scattered photons are distributed uni- 
formly over the whole 47r solid angle. Note that this is 



in contrast to the weak pumping case < 7, in which 
the first-order correlation function does show directional 
behavior as 



1 



+N{N -l)\{S+)f cos[(fc-fcL)f,,]} , (9) 



where 



(5+) = 



7' + A' 



2(72 + A2 + 172) ' 
z»(72 + A2) 



(10a) 
(10b) 



D. Second order correlations 

Our main observable in the following will be the the 
normalized second-order correlation functions g^'^^{T), 
which is measured by two detectors positioned at i?i and 
R2. It is defined as 



(2)/ D 3\ Gmn{T,Ri,R2) 

g^4(T, i?i,i?2) = : 

(i?l)/„(i?2) 



(11) 



i.e., as the unnormalized second order correlation func- 
tion 

G^^iiT,R,,R2)={al{Ri)aUT,R2) 

X o„(r,^2)am(^i)) , (12) 

normalized to the intensities Im{Ri) and /„(i?2) defined 
in Eq. ([6]). The quantity gmn for {m, n} € {C, R, L} 
can be interpreted as a measure for the probability of 
detecting one photon emitted in mode m and another 
photon emitted in mode n with time-delay r. Particu- 
larly, g^'^^T = 0) describes the photon statistics (e.g., 
sub/super-Poissonian), whereas g'^^^T ^ 0) indicates 
photon bunching or antibunching. 

To calculate the correlation function, we assume laser 
driving on resonance (9 = 7''/4), and a large atomic en- 
sembles {N ^ 1) such that the secular approximation is 
valid. We also assume that all possible pairs of atoms 
contribute equally to the second-order correlation func- 
tions. This assumption is valid as long as the angle be- 
tween the wave vectors of the incident laser and the scat- 
tered photons is small, that is {(j), cpo} should be of order 
of few degrees (see Fig. [T|) . Finally, we for the moment 
consider a single interparticle distance vector fji for all 
pairs only, but this restriction will be relaxed later on. 
Based on these assumptions, we found that the correla- 
tion and cross-correlation functions of photons scattered 
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into the different spectral bands can be represented as 

g^cU^, Ri,R2) = 1 + 2 cos(<5i) cos(52)e-'^" , 

= 1 + [cos(<5+) + cos((5_)] , (13a) 

= l + cos((5_)e-37r^ (13^) 

= 1 + cos((5+)e^3'^^ , (13c) 

g^^U^,Ri,R2)^g%{T,R,,R2) 

= 1 forXe{L,i?}. (13d) 

Here, 



S+ = Si + 62 , 
S- ^ 5i - $2 , 



(14a) 
(14b) 
(14c) 



with kg being the wave- vector of the photon s scattered in 
direction Rs (s € {1,2}). Note that these results differ 
from the corresponding expressions for a single pair of 
atoms or a regular structure of atoms [2^ . The reason 
is that for a single pair, one chooses iV = 2, whereas in 
deriving Eqs. p3)) the assumption 3> 1 was used such 
that terms proportional to 1/A^ can be neglected. 



2. Averaging over a spherical volume 

This method is a modification of the first method in 
which we change the averaging over the radial coordinate 
to range from to twice the volume radius R: 



(...) 



1 



conf 



8ttR 



2B. 



(16) 



This averaging procedure marks the opposite extreme of 
the averaging in Sec. HIE II in that a constant distance 
distribution over all possible distances is assumed. In 
the numerical examples below, we use R = IOOAl- 



3. Numerical sampling 

The first two averaging methods made assumptions on 
the distribution of distances between the atoms, which 
may not be fulfilled in a gas of atoms. To investigate the 
effect of this further, in a third averaging procedure, we 
randomly place a number of atoms TV in a cubic volume 
of side length 2R and then calculate the quantities in 
Eqs. (jl3p for this sample of atoms. In the numerical 
examples below, we use R = IOOAl and N = 300. 



E. Ensemble averaging 

To estimate the signal obtained from a cloud of ran- 
domly distributed particles, the correlation functions 
gxyi'T: RitRt) have to be averaged over the different in- 
teratomic distance vectors Vij in the cloud. Since in gen- 
eral it can be expected that the averaging procedure af- 
fects the result, we compare three averaging procedures. 

1. Averaging over a spherical shell 

The first averaging procedure we consider has been 
suggested in previous works on scattering from dilute 
gases, see, e.g., It consists of (i) an isotropic av- 

erage over the relative orientation n of the atoms over 
the unit sphere, followed by (ii) an average of the inter- 
atomic distance Tji = \ over an interval of order of the 
laser wave-length, around their typical distance I : 

{■ ■ ■ )conf = .0 / drji / diln ■■■ . (15) 

(,47rj Jl-2Tr/kL J 

It should be noted that the restriction to the radial av- 
eraging over a range of order of the wave length is not a 
priori justified; it is however motivated by the fact that 
the distance distribution of particles in a gas is peaked 
around a mean value, even though the distribution gen- 
erally is broad. In the numerical examples below, we use 
I = 20Xl. 



III. RESULTS 

As a first result, we note from Eqs. ([T^ that the sec- 
ond order correlation function can exceed unity in sev- 
eral cases together with bunching, which means that 
there is an enhanced probability to generate photons in 
pairs. For example, the second-order correlation func- 
tion for the central-band photons g^^ has maxima when 
ki + k2 — "^kL [see Eq. IT^a)]. i.e., in forward direction. 
Similar results are observed for the sideband spectral 
lines. For example, pairs of photons can be generated 
in which one photon is emitted in the left spectral band, 
whereas the other is emitted from the right spectral side- 
band, see Eq. lT^ c). In contrast, photon pairs originating 
from the same sideband (left or right) are most likely for 
ki = k2 independent of the direction of the wave vector 
and thus do not show directional behavior, see Eq.fT^Tb). 
Finally, there is no enhanced probability for detecting 
pairs of photons in which one photon is generated in the 
central band and one in either of the sidebands, and the 
corresponding correlation function does not exhibit any 
directionality, see Eg. [T^d). 

Regarding the intensity, the unnormalized second- 
order correlation functions are proportional to the num- 
ber of particles in the ensemble squared, i.e. G^'^^{0) oc 
N\ 

The obtained results for the normalized second-order 
correlation function 17^^(0) for photons emitted in the 
central band after the configuration average are shown in 
Fig. [21 It can be seen that even after the averaging, the 
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FIG. 2: (Color online) Normalized and configuration aver- 
aged second-order correlation function g^ci^) between two 
photons emitted from the central spectral band. In (a), the 
correlation function is shown for pairs of photons emitted in 
the same direction (0o = 0, see Fig.[l|, and plotted as a func- 
tion of the emission direction (j>. In (b), the correlation func- 
tion is plotted as a function of the opening angle 4>o between 
the two photons, which corresponds to the case of two distinct 
detectors. The two photons are measured at positions sym- 
metric with respect to the incident laser field direction, i.e., 
<j!> = 0. The blue solid line shows averaging over a spherical 
shell, the dashed red line averaging over a spherical volume, 
and the green dotted line the numerical sampling. 



correlation function exhibits a sharp peak around direc- 
tion of the the incident laser wave- vector fc^, indicating 
super-bunching. Fig. [H^a) shows the case with a single 
two-photon detector (00 = 0) for different emission direc- 
tions (t>. Figure[2jb) depicts the same correlation function 
but with two single-photon detectors placed symmetri- 
cally with respect to the laser wave-vector direction 
given hy (f) ~ {). Note that the correlation functions for 
two photons emitted both from the left or both from the 
right spectral sideband do not show a directionality in 
space. 

Fig.[3]shows the corresponding results for the two pho- 
ton cross correlation with one photon emitted in the left, 
and one in the right spectral sideband. Again, we find 
maxima at i^o = for two individual detectors placed 
symmetrically around the incident laser direction. But 
in contrast to the central band correlation function, this 
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FIG. 3: (Color online) Same as Fig. (2] except that the nor- 
malized and configuration averaged second-order correlation 
functions ff^^^^^^j (0) between one photon emitted in the left 
and one photon emitted in the right sideband is shown. 



maximum is not peaked, but rather broad. If a two- 
photon detector is used, a narrow maximum is observed 
in the forward direction, see Fig. [H^a). Interestingly, in 
this case, depending on the precise positioning around 
the forward direction and on the averaging procedure, 
also sub-Poissonian photon-statistics can occur. 
We further define the Cauchy-Schwarz parameter 



5il(0) 



(17) 



which relates the correlation between photons emitted 
into individual modes to the cross-correlation between 



photons emitted into two different modes j3j|. If x < 1 
the Cauchy-Schwarz inequality is violated. Figure |3] 
shows violation of the Cauchy-Schwarz inequality for 
photons scattered into the side-bands. 



IV. DISCUSSION AND SUMMARY 

An interpretation of the results obtained in the pre- 
vious section can be found in the semiclassical dressed- 
state picture. Suppose one laser photon is absorbed and 
the atom is in the dressed-state |2). Then one option 
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FIG. 4: (Color online) The Cauchy-Schwarz inequality x 
characterizing the cross-correlations of photon pairs with one 
photon emitted from each sideband. The photons are de- 
tected symmetrically around the incident laser field direction 
(0 = 0), and the result is shown as a function of the angle 
between the wave vectors of the two emitted photons cj>o. The 
Cauchy-Schwarz inequality is violated for x < 1. The curves 
are as in Fig. [5] 



will be a spontaneously emitted photon on the |2) |2) 
transition followed by a further absorption of a laser pho- 
ton during time- interval 17^^ and subsequently decay on 
the same dressed-state transition. Similar effect occurs 
on the dressed-state transition |1) o |1) when initially 
the emitter is in the |1) dressed-state, and both pro- 
cesses contribute to the central-band scattering where 
ki + k2 2kL. Another option will be a spontaneously 
emitted photon on the |2) |1) dressed-state transi- 
tion followed by an absorption of a laser photon on the 
|1) —5- |1) transition and further spontaneously decay on 
the |1) 1 2) dressed-state transition. Together with the 
reverse process where a spontaneous emission occurs on 
the |1) — i> 1 2) transition (if initially the particle is in the 
|1) dressed-state) followed by a laser absorption on the 
|2) 1 2) transition and a subsequently decay on the 
|2) — J> |i) dressed-state transition, these effects describe 
photon correlations among the right and left side-bands, 
or vice versa and, again, fci -f fc2 ~ 2A;l, and therefore, 
we have a directional-dependent photon distribution. In 



contrast, all other scattering processes, i.e., two-photon 
spontaneously emission on the same side-band or cross- 
correlations involving a side-band photon and a central 
one do not fulfill the relation ki + k2 ~ 2kL, that is 
^1 + ^2 ^ 2fcL, and thus, in these processes the photon 
distribution will be uniformly in space after a configura- 
tion averaging. 

Comparing the results for the three averaging proce- 
dures in Figs.[2]|ll it can be seen that the averaging strat- 
egy can substantially affect the result quantitatively, even 
though the qualitative features of the obtained results re- 
main similar. In particular the averaging over a spherical 
volume and the numerical sampling in a cubic volume of 
comparable size agree reasonably well. A notable dif- 
ference between the three methods is the sub-Poissonian 
statistics found in Fig. [3l^a), which is pronounced only in 
the case of averaging over the spherical shell. The reason 
is that the fringe pattern seen when plotting against (fi 
washes out if the correlation function is averaged over a 
larger range of distances. 

In summary, we discussed the second order correla- 
tion function for light emitted from a strongly and near- 
resonantly driven dilute cloud of atoms. Because of the 
strong driving, the fluorescence spectrum separates into 
distinct peaks, for which the spectral properties can be 
defined individually. We have shown that the second- 
order correlations for various combinations of photons 
from different spectral lines exhibit bunching together 
with super- or sub-Poissonian photon statistics, tunable 
by the choice of the detector positions. Furthermore, 
we demonstrated that the Cauchy-Schwarz inequality is 
violated for photons emitted from the two spectral side 
bands. The emitted light intensity is proportional to the 
square of the number of particles, and thus can poten- 
tially be intense. Interesting applications might arise if 
the presented method is applied to generate correlated 
photon pairs in the x-ray domain [ssl - 37 1 . 
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